Optimal correlations in many-body quantum systems 
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Information and correlations in a quantum system are closely related through the process of mea- 
surement. We explore such relation in a many-body quantum setting, effectively bridging between 
quantum metrology and condensed matter physics. To this aim we adopt the information-theory 
view of correlations, and study the amount of correlations after certain classes of Positive-Operator- 
Valued Measurements are locally performed. As many-body system we consider a one-dimensional 
array of interacting two-level systems (a spin chain) at zero temperature, where quantum effects are 
most pronounced. We demonstrate how the optimal strategy to extract the correlations depends on 
the quantum phase through a subtle interplay between local interactions and coherence. 
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The relation between correlations and measurement 
is well known in quantum metrology, where the opti- 
mal measurement strategy to extract information has 
been thoroughly investigated PQ [2]. In that context 
Positive-Operator- Valued Measurements (POVMs) and 
Informationally-Complete (IC) measurements are of par- 
ticular importance since, contrary to simple projective 
measurements, they allow a complete tomography of the 
quantum state [5] . They have been playing an important 
role also in foundational aspects of quantum mechanics, 
quantum information science, as well as in the physics of 
dissipative systems [IHZI- 

In this Letter we aim at establishing a bridge between 
quantum metrology and quantum many particles physics. 
We consider subsystems A and B in a many-body ground 
state, and analyze the correlations resulting from POVMs 
performed on one of them, say B. Emphasis will be de- 
voted to the optimal correlations, namely the maximal 
amount of correlations established between A and B. We 
observe that performing a POVM on a given physical sys- 
tem is equivalent to performing a projective measurement 
on an enlarged system where the original one is coupled 
with a given "ancilla" (Naimark's theorem). Such an an- 
cilla can be an appropriate subsystem, and then analyz- 
ing correlations induced by a POVM on a local degree of 
freedom, say B, is an effective way to study correlations 
of higher order (spin-block correlations). Equivalently, 
the ancilla can be a suitable environment entangled to 
the system, and then correlations can give precious in- 
formation on the decoherence of the state of the local 
constituent A. 

The total amount of correlations in any bipartite 
(mixed) quantum state pab is given by the mutual in- 
formation: Tab = S{pa) + S{pb) — S{pab), where 
S{p) — — Tr[plog2/5] is the von Neumann entropy. A 
central quantity we will address below is the quantum 
conditional entropy Sc , quantifying the ignorance about 



the composite system AB, once subsystem B has been 
measured with a generic POVM {-Bfc}: 

Sc{pab\{B,})^J2p'^^(PaI)- (1) 

k 

Here denotes the state of the composite system AB, 
conditioned to a given outcome of Bk- Pab ~ p^t'-^ ® 
Bk) Pab] with I denoting the identity operator on the 
subsystem A and pk = Tr[(I (g) Bk) pab]- The global 
amount Cab of optimal (classical) correlations between 
constituents A and B is established after applying a set 
of measurements on B that disturbs the least the part A: 

Cab - max [5(pa) - ScipAB\{Bk})] , (2) 

where the maximization is with respect to the mea- 
surement strategy. The quantum discord [5] quantifies 
the amount of quantum correlations and is defined as 
the difference between total and classical ones: Qab = 
Iab — Cab ■ The maximization in Eq. ([2| is generally a 
daunting task, since the optimization procedure has to 
be performed on the whole set of possible POVMs. 

We apply the above notions to the case where A and 
B are individual spins of a quantum spin chain, and con- 
sider both von Neumann projective measurements Mproj, 
and generalized POVMs Mpovm We design a strat- 
egy to exploit the information input given by the physical 
system hosting the two spins. Namely, we assume that 
the symmetry of the POVM is fixed by the symmetry of 
the local interactions occurring in the physical system. 
However, we shall see this is not enough to optimize cor- 
relations, as the coherence of the many-body system is 
going to play a major role. 

Models and measurements. — As many body system 
we consider a one dimensional array of spins 1/2 inter- 
acting anisotropically along the three spatial directions 
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with interaction strengths Jx, Jy, Jz, and subjected to a 
uniform external field h. The Hamiltonian reads 



Jy 



(3) 

where af {a = x,y,z) are the Pauli matrices on site 
i. Hereafter we set \Jx\ = 1 as the energy scale and 
work in units of /i = 1. At zero temperature different 
quantum phases can exist, separated by Quantum Phase 
Transitions (QPTs) [10]. Moreover a completely factor- 
ized ground state may occur at a specific value of the 
field h f [TT] . For xyz spin chains of Eq. Q this is given 
by: hf = 2^{l-i){Jy-J,) m- 

We will discuss Eq. ([s]) in the following cases: I) The 
ferromagnetic Ising chain {Jy — Jz = 0, Jx — —1), which 
undergoes a QPT at he — 1, and factorization aX hf — 0. 
It can be experimentally realized with the magnetic com- 
pound Co Nb2 Oe [13] . II) The non-integrable antiferro- 
magnetic xyx model { Jx = Jz = l)j with Jy = 1/4 (this 
is the case experimentally realized with Cs2 Co CI4 [14j). 
Such model displays a QPT at he — 3.21 and factoriza- 
tion aX hf c^i 3.16. Ill) The antiferromagnetic anisotropic 
XX z Heisenberg chain {Jx = Jy = 1) with = 1/2. At 
zero field it presents a critical xy phase with quasi-long 
range order (quasi- Iro) for \Jz\ < 1; this is separated by 
two classical phases with QPTs at Jz — ±1. For h ^ 
the xy phase is a strip in the phase diagram, eventu- 
ally turning into polarized phases for sufficiently strong 
magnetic field (the factorization phenomenon degener- 
ates in the saturation occurring as a first order tran- 
sition). Despite local interactions are clearly different, 
both the quantum Ising and xyx models display an Ising- 
like QPT with Z2-symmetry breaking; the xxz model, in- 
stead, is characterized by a critical phase without order 
parameter. 

We first deal with standard projective measurements 
Mproj = {B±} along the field direction, defined by 
B± — 5(1 ± (T^). Then we engineer a more sophisticated 
set of POVMs, such that the symmetry of the measure- 
ment keeps track of local interactions between the spins. 
Specifically, we look at the interactions Jx, Jy, Jz enter- 
ing Eq. (|3|, and design the following Mpovm — {Bk}'- 
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A; = 1...4, 



(4) 



and Ofc is such that 
a.{Jx,-Jy,-Jz), as = 



where a = (ct^, , 

Si = a{Jx,Jy,Jz), 0,2 : 

a{-Jx,Jy,-Jz), 04 = a{-Jx,-Jy,Jz) aud = 
Jx + Jy + Jz (see Fig. jlj). For generic Ja, Mpovm will 
be denoted as Coupling-oriented Informationally Com- 
plete (CIC) POVM. The choice of the vectors dk in 
Eq. Q refiects the symmetry of the Hamiltonian by 
changing J^ - Ja by UUW, with [>" = U^^2^+l 



J. 



In the isotropic case Jx — Jy — Jz, the POVM in Eq. H\ 



FIG. 1: (color online). The four vectors entering the POVM 
measurement of Eq. Q. They point from the center to non- 
adjacent corners of a parallelepiped with edges fixed by the 
anisotropic interaction occurring into the system. 



degenerates in a Symmetric Informationally Complete 
(SIC) POVM [5]. We comment that all the Bk do not 
depend on the external field explicitly. We shall see that 
such a dependence enters in a subtle way related to the 
macroscopic order of the system. 

In order to compute the amount of correlations be- 
tween any two spins at distance r = \A — B\, one needs to 
access the single- and two-spin reduced density matrices 
Pa and pab- Hereafter we focus on the symmetry-broken 
ground states of the Hamiltonian in Eq. ([3|, which is 
symmetric under a global phase fiip along the z-axis [15j . 
We observe that the generic two-site reduced density ma- 
trices of such states is beyond the so called "X-state" 
structure (emerging in symmetric states, and for which 
expressions for quantum and classical correlations are 
known [16 ). Therefore, in the present case, in principle 
the optimal correlations might be achieved beyond pro- 
jective measurements [T7]. To access all the required two- 
point correlators, we resort to the density matrix renor- 
malization group approach with open boundary condi- 
tions [TH]. We consider sufficiently long chains, such to 
reduce unwanted edge effects and to approach the ideal 
thermodynamic limit. For the Ising and the xyx model, 
we add a small longitudinal field hx ~ 10~^ along the xy 
plane, to ensure the Z2-symmetry breaking. 

Comparison between different measurement 
strategies. — We start our analysis by presenting 
results obtained for the quantum conditional entropy 
S^^ in Eq. ([T]), probing how the local interactions affect 
the measurement, without any further optimization. 
Figure [2] displays Sq'^ for two neighboring spins re- 
spectively for the Ising model, the xyx model and the 
xxz model in a transverse field h. The total amount of 
correlations reflects the main properties of the ground 
state: in particular the peaks denote QPT points that 
are associated to a divergence in their first derivative, 
while factorization fields are marked by the vanishing 
of all correlations. In all the three considered spin 
systems, measurement performances decrease from the 
CIC POVM to the SIC POVM and to the projective 
measurement along the computational basis z. 

Much larger amounts of correlations can be achieved by 
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FIG. 2; (color online). Gonditional entropy 5^^' for the three 
cases: I) Ising, II) xyx and III) xxz model in an external 
field [see Eq. (|3|]. Correlations are considered between two 
nearest neighbors at the center of the chain (r = 1); the length 
of the chain is L = 200. The various curves correspond to 
different measurement strategies on subsystem B: standard 
projective measurement along the z direction (black circles), 
SIC POVM (red squares), as well as IC POVM set by the 
specific interactions of the model (blue diamonds). 



FIG. 3: (color online). Same as in Fig. [2j but for the clas- 
sical correlations C'^' optimized over the given measurement 
strategy. In the first three cases, projective measurement, SIC 
POVM and IC POVM (see the symbol pattern of Fig.^, the 
optimization in Eq. ^ is performed on the rotation angles 
{0,(j)) of the Bloch sphere of subsystem B. The green stars 
refer to IC POVM optimized by varying the three parameters 
Ja defining the direction of the vectors a^- 



performing suitable optimization strategies on the mea- 
surements considered above. In the following we apply 
two different kinds of optimization: i) We rotate the di- 
rection of the elements Bk of the projective measurement 
Mproj and of the POVM Mpovm on the Bloch sphere, by 
keeping the angles between the vectors Oj constant (this 
corresponds to a rigid rotation of the experimental ap- 
paratus) ; optimal correlations are thus obtained by max- 
imizing over the angles (0, 0) entering the rotation, ii) 
In the case of CIC POVM, we independently vary the 
three parameters Jx,Jy,Jz, defining the direction of the 
vectors dk in Mpovm (see Fig.[l]) [T5]. 

In Fig. [3] we display the classical correlations between 
two neighboring spins for the three considered models, by 
adopting the optimizations discussed above. Similarly to 
the quantum conditional entropy, C'""^^^ displays a no- 
ticeable dependence on the model and on the magnetic 
field. While in the Ising model the CIC angle-dependent 
POVM and projective measurement give the same an- 
swer, for the xyx and xxz model the angle-dependent 
strategy is not optimal and is outperformed by the pro- 
jective measurements. The 3-parameters POVM opti- 
mizations provide the same correlations of the projec- 
tive measurements in the disordered regions and in xxz 
model; however, where the order parameter (cr^) 7^ 0, 
they are still worse than the projective measurement. 

A similar analysis of correlations for r > 1 strongly 
suggests that the effect of different measurement strate- 
gies at long ranges vanishes everywhere but close to the 
quantum critical points, where the correlation functions 
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FIG. 4: (color online). Quadratic difference between the op- 
timal correlations obtained following the four different mea- 
surement strategies M described above (see also Fig. [sjl: 
= YIm (C'''' - {C^''^)f- Here (C'''') denotes the average 
correlation with respect to the four M. Data are shown for 
the xyx model as a function h and for different r = |A — _B|. 



decays algebraically with r (Fig. |4]). 

In the disordered phase of the Ising and the xyx model, 
where the order parameter (ct^) = 0, as well as in the xxz 
model, Oopt and 4>opt are fixed to a value independent by 
h (see Table |l]) . By analyzing the rotated measurements 
Bk{0,(j>) it turns out that projective measurements can 
achieve a local measurement along the eigenvectors \v) 
of (Jioc = Jx<y^ + Jy<^^ + Jz^^ , fixed by the system 
Hamiltonian. For both the xyx and Ising model, opti- 



4 



mal correlations are attained by using projective mea- 
surements along the x axis: B±{-k/2,Q) = a^). 
This reflects the 2'2-symmetry cr^ of the para- 

magnetic phase. On the other hand, the operators of 
each optimized 4-elements POVM can be written as: 
Bk{6opt,(t>opt) = where is of the type 

W) - a^/r^e^'-^-'i t) + \/rT^| ;), with a e c. 

It turns out that IV'fc^*) = with (l>opt = and 
S,k — 'C'"'' = g«arg(a^-(-ia„) projective measurements. 
For CIC measurements Cfc = C^''^ but (f) ^ 0; for SIC- 
POVM ^fc ^ ^("l We note that for large h, where the 
state is nearly fully polarized along z, correlations are 
vanishing and therefore measurements along any direc- 
tion are optimal. The CIC with three-parameter opti- 
mization leads to optimal correlations in the disordered 
phase. 



In the symmetry-broken phase, optimal angles 
(^opt, 0opt) for both Ising and xyx models display a non 
trivial dependence on the order parameter, as visible in 
Fig. [5] Except for a region close to the QPT, we fitted 
our results using the formula 9 opt — B — [a^)^ + k, 
where the parameters A, B, k are model-dependent, while 
we imposed n = 8, inspired by the linear variation of 9opt 
with h that we observed in the Ising model, and by the 
characteristic exponent = 1/8 of its order parameter 
at the QPT (see the upper inset of Fig. [s]). In some 
cases we found significant deviations from such depen- 
dence (see Supplementary Material for a detailed discus- 
sion of the optimal parameters) . Peculiar behaviors arise 
close to the QPT, where dramatic changes appear, and 
to the factorization points, where an extremal value is 
reached. Such behavior is consistent with the interpre- 
tation of the factorization phenomenon as a "correlation 
transition" resulting from a competition between parallel 
and antiparallel entanglement [20l [21] : Optimal correla- 
tions arise from the balance of two optimizations involv- 
ing the parallel and antiparallel entangled components 
(both present in the ground state). The two entangle- 
ment components switch each other and an extremal op- 
timal angle is reached at hf. 

It is interesting to compare the optimal angles in the 
Z2 symmetry-broken phase with those in the xy gapless 
phase of the xxz model, where the order parameter is 
vanishing in a non trivial way because correlations de- 
cay algebraically. For both projective and CIC measure- 
ments, the optimization in such phase is characterized by 
a fixed value of 9opt, '^4'opt, thus indicating that, because 
of quasi-long range order, any preferential measurement 
direction is not unique in the phase. Such scenario is con- 
firmed by the three-parameters optimized CIC POVM 
(last row of Table |l]). 

Discussion. — We analyzed spin-spin correlations that 
are established after performing a local measurement on 
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FIG. 5: (color online). Optimization parameters for the 
nearest-neighbor correlations in the xyx model, after a pro- 
jective measurement. The fit in the upper inset (dashed blue 
line) is Qopt = 0.824-0. 709\/0.0769 - x^. The other optimiza- 
tion strategies are discussed in the Supplementary Material. 
In the lower inset we display the order parameter (a^) as a 
function of h, as extracted from numerical simulations. 



one of the two spins in the ground state of a quantum 
spin chain. We considered projective measurements as 
well as symmetric IC POVMs; furthermore we engineered 
coupling-oriented IC POVMs with the aim to shed light 
on how the optimal measurement can be performed a 
priori, once a certain knowledge on the system has been 
previously acquired. The measurement protocols were 
first tested regardless to adjustable parameters, by look- 
ing at the conditional entropy. Then we focused on the 
possibility to adjust the measurement on the basis of lo- 
cal interactions. Interactions, coherence and symmetry 
fix the "optimal flow" of information; the optimal strat- 
egy to extract the correlations eventually depends on the 
quantum phase. 

Specifically, an analysis of the quantum Ising and the 
xyx model, both sharing the same kind of Z2-symmetry 
breaking QPT (even if with very different local interac- 
tions), showed that in the ordered phase the optimal cor- 
relation follows the global order, in the sense that the pa- 
rameters characterizing the optimal measurement strat- 
egy vary with the exponent /3 of the order parameter. 
Such a result could be useful to hcuristically dictate the 
optimal measurement strategy also for higher order cor- 
relations (between spin-blocks), where the optimization 
protocol is not practicable. On the contrary, in the dis- 
ordered phase local interactions fix the optimal strategy, 
in the sense that optimal correlations are attained by ful- 
filling a local requirement of projecting along cr;oc- The 
results on the xxz model further support this scenario: 
optimal correlations are obtained for measurements re- 
specting the in-plane symmetry of the model, for any 
fixed direction in the xy plane (there is no preferential 
direction because of quasi-lro in the xy critical phase) . 

Finally we analyzed correlations at long ranges show- 
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Ising, h > he 


xyx^ h > he 


XX z, quasi- Iro 


Proj. 


01 = 7r/2; = 

02 = 0i; </>2 = 


01 = 7t/2; <t,i=0 

02 = 01, (t>2 = TT 


= TT/2;\f(t) 


CIC 


01 = 0; V(/.i 
02 — it; (j)2 — 0,n 
Wa; (jj-i = TT /2 + kn 


0k = tt; 
(pk ~ 1-478 + kTT/2 


e = 0; V0 


SIC 


0k = tt; 
cj>k « 0.393 + kTT/2 


0k = tt; 
4>k « 1.152 -|-fc7r/2 


0k^ ~ 0.955; (j}k^ = 37r/4 + km 
0k2 ~ 2.185; (pk2 ~ 3.92 ferr 


CIC 3-par 


Jx — I5 Jy — Jz — 


tlx — 1) Jy — Jz — 


= 0; V(Ji:, Jy) 



TABLE I: Optimization angles 0opt G [0, vr] and (popt G [0, 2tt) in the disordered phases with vanishing order parameter, k is 
any integer positive number. In the last line: optimization directions Ja for CIC measures. 



ing that, near the QPT, long-range correlations are 
strongly affected by the measurement strategies (see 
Fig. [4]). In the gapped phase any measurement strat- 
egy produces the same result, on a length-scale where 
the correlation functions themselves are sensible. 

Given the relation between optimal correlations and 
quantum discord, our results could be important in 
many-body implementations of quantum information 
protocols. We also comment that, being a single-spin 
state fully accessible through Eq. Q, our scheme pro- 
vides an effective strategy to perform state tomography 
of one of the two spins (a notoriously challenging problem 
in quantum magnets). 
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Supplementary material. 

Here we report a discussion on the optimal angles 
{0opti 4>opt), as well as on the optimal vectors {Jxj Jy, Jz)i 
which maximize the correlations in the symmetry-broken 
ordered phase. Namely, we focus on h < he in the Ising 
and the xyx model. Only the case r = 1 of two nearest- 
neighbor spins is addressed. 



ISING MODEL 



the figure (up to numerical precision). Here we select 
this point for continuity with the {9i = tt; 0i 0.393) 
optimal point in the disordered phase [see Table (1) in 
the main text]. 

For C-IC optimized on the three directions {Jx, Jy,Jz), 
the optimal point is found to be at Jy = 0, while Jx and 
Jz vary with h according to Fig. [6]:. Away from the 
critical point, we can fit Jz linearly with h. Note that 
the value of Jx is fixed by the normalization constraint 
llflfcll — 1, and its variation is much smaller on the scale 

of Jz. 



For both projective and C-IC measurements there are 
two optimal couples of angles. The first pair is deter- 



mined by (j)i^Jt = and 9),p[{h = 0) = 7r/2, with 9),p't oc h 
(Fig. |6^), as long as h increases until the critical point 
he — I (note however the dramatic changes in 9opt close 
to he). The other pair of optimal angles is (tt — 9*^^^; n). 
The approximately linear behavior of 9opt with respect 
to the field h results in a fitting law analogous to Eq. (6) 
in the main text, namely: 



3(1)/ 



3(1) 



^opt 



= AyjB - ((7^)" + fc . 



In the specific we found A w 0.34 and i? = 1; this is 
because, for the Ising model, (ct^) = \\ — ]r?\^l^ . 

For SIC POVM 9opt displays corrections to the linear 



behavior with h, while 



'->opt 



const., as compared to 



the scale of variation of 9opt (Fig. [6j3). We point out 
that there are other optimal points which maximize the 
correlations to the same amount of the ones depicted in 



XYX MODEL 

As in the Ising model, for the projective measures we 
found two optimizing points. The dependence of 9opt 
with h for one of the two points is depicted in Fig. (5) of 
the main text, while the corresponding axial angle 0opt = 
0, independently of h. The other pair of optimal angles is 



(5) ('"' ~ ^opt' '"')' ill complete analogy with the Ising model. 
We fit 9opt as a function of the order parameter (cr^) 
as in Eq. ([5 1; the fitting parameters are A = 0.709 and 
B = 0.0769 see dashed blue curve in Fig. (5) of the main 
text] . 

For the other measurement strategies we observed de- 
viations from this behavior and we could not operate 
such a fit. The dependence of the optimal angles on h is 
displayed in the various panels of Fig. [Tj 
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FIG. 6: Optimization parameters for the nearest neighbor correlations in the Ising model. The different panels correspond 
to different optimization strategies, as discussed in the main text. Red dotted-dashed lines are linear fits of data, while blue 
dashed ones are power-law fits; fitting parameters are indicated in the respective legends. The tilde in panel a) stands for a 
rescaling of Oopt such that 6opt{h = 0) = 0. 
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FIG. 7: Optimization parameters for the nearest neighbor correlations in the xyx model. The different panels correspond to 
different optimization strategies, as discussed in the main text. The insets display how the optimal angles are pronounced 
non-linear functions of (cr^). 



